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High-accuracy scaling exponents in the local potential approximation 
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We test equivalences between different realisations of Wilson's renormalisation group by comput- 
ing the leading, subleading, and anti-symmetric corrections-to-scaling exponents, and the full fixed 
point potential for the Ising universality class to leading order in a derivative expansion. We discuss 
our methods with a special emphasis on accuracy and reliability. We establish numerical equivalence 
of Wilson- Polchinski flows and optimised renormalisation group flows with an unprecedented accu- 
racy in the scaling exponents. Our results are contrasted with high-accuracy findings from Dyson's 
hierarchical model, where a tiny but systematic difference in all scaling exponents is established. 
Further applications for our numerical methods are briefly indicated. 



I. INTRODUCTION 

Renormalisation group (RG) methods provide 
powerful tools in the computation of universal 
scaling exponents or anomalous dimensions in 
quantum field theories Modern functional RG 
approaches are based on Wilson's idea to integrate- 
out momentum modes within a path-integral 
representation of the theory (2j. A particular 
strength of these methods is their flexibility when 
it comes to approximations. Wilsonian flows can be 
implemented in many different ways, which helps 
adapting the technique to the problem at hand @]. 

A useful systematic expansion scheme is the 
derivative expansion, which, to leading order, is 
known as the local potential approximation (LPA). 
For the study of critical phenomena, the derivative 
expansion is expected to provide good results as 
long as quantum corrections to propagators and 
anomalous dimensions remain small [4|. Important 
results in LPA have been accumulated over the 
years based on different implementations of the 
RG including the Wilson-Polchinski RG @, i|, the 
functional RG [2] for the effective average action 
and optimised versions thereof [1], and Dyson's 
hierarchical RG 0, [l(| ■ Interestingly, and despite 
qualitative differences in the respective formalisms, 
these approaches are closely related in LPA. In 
fact, it has been established by Felder [Tl[ that the 
Wilson-Polchinski RG and the hierarchical RG are 
equivalent. More recently, following a conjecture 
first stated in 0, [H[ , the equivalence between the 
Wilson-Polchinski flow and an optimised functional 
RG has been proven by Morris [14j. Consequently, 
universal scaling exponents from either of the three 
approaches should be identical. 

Quantitatively, these equivalences have been put 
to test for 3d critical scalar theories [13, UM . While 



scaling exponents from either approach agree at the 
order 10 -4 , an unexpected disagreement starting 
at the order 10 -5 was found. Presently, however, 
this last observation solely relies on results from 
the hierarchical RG and the optimised RG, where 
exponents have been obtained with sufficiently 
high accuracy. Therefore it becomes mandatory to 
obtain results from the Wilson-Polchinski RG with 
a higher precision, in order to confirm or refute the 
above-mentioned discrepancy. 

In this paper, we close this gap in the literature 
and study both the Wilson-Polchinski RG and the 
optimised RG for the 3d Ising universality class with 
an unprecedented accuracy. To that end, we in- 
troduce several methods to solve non-linear eigen- 
value problems with special emphasis on the nu- 
merical reliability. We expect that these techniques 
are equally useful for other non-linear problems in 
mathematical physics. The outline of the paper is as 
follows. We introduce the relevant differential equa- 
tions (Sec. [Til and discuss our numerical methods 
(Sec. IIII|) and the error control (Sec. IIV[) . Results 
for the fixed point solution and scaling exponents 
are discussed (Sec. |Vj and compared with the hier- 
archical model (Sec. lVI|) . We close with a discussion 



and conclusions (Sec. IVII 



II. LOCAL POTENTIAL APPROXIMATION 

We restrict ourselves to a real scalar field theory 
in three dimensions - the Ising universality class - 
and to the leading order in the derivative expansion, 
meaning that the scale-dependent effective action 
is approximated by Tfc = J d 3 x[^d^(j)d^(l) + I4(</>)]. 
Here, k denotes the scale of the Wilsonian mo- 
mentum cutoff. We provide the flows in terms of 
the dimensionless potentials u(p) — Vk(<fi)/k 3 with 
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p = ^<j> 2 /k, and v(<p) = Vk((f>)/k 3 with tp = 6/Vk. 
Within an effective average action approach [Tj , the 
flow equation for u(p) depends on the infrared mo- 
mentum cutoff. For an optimised choice of the latter 
Q, and neglecting an irrelevant (field-independent) 
term, it is given by 

d t u = -3u + pu + — — ^— -. (1) 

1 + u + 2pu" 

Here, t = ln/c denotes the logarithmic scale param- 
eter. At a fixed point dtv! — 0, the potential obeys 

du" 
dp = 



(1 



2p u' 



3u" 
~2p~ 



(2) 



There exists a unique non-trivial and well-defined 
(finite, no poles) solution to (2} which extends over 
all fields. In terms of the potential v(tp), the corre- 
sponding equations are 



d t v 



-6v + 2 t fi v 



1 



(3) 



(4) 



1 + v" 

^ = \{<pv"-hv>){l+v"f 

For the Wilson-Polchinski RG, we follow an analo- 
gous ansatz for the effective action. In terms of the 
potential u(p) and with t = In k wc find 

dtu = —3u + pu—u' — 2pu" J r2p(u') 2 (5) 

du" u" u' {u'f n , „ 3u" , . 
— = — -- + ^-L + 2u'u" - — . (6) 
ap 2 p p 2p 

The corresponding equations for v((p) are 

1 



-3v+ ^ipv' -v" + {v') 2 , (7) 



[tp v" - 5v') + 2v' v h 



(8) 



d t v 

dv" 1 
d^" ~ 2 

Numerical solutions for the fixed point potentials 
and their derivatives @, (UJ) and ([5]) are displayed 
in Figs. [[J H El El and below. 

The potential v(ip) in ([3]) is related to the Wilson- 
Polchinski potential i>wp (vwp) in J7J) by a Legendre 
transform [141 ] . 
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WP ! 



(9) 
(10) 

up to an irrelevant constant (the absolute value 
of the potentials at vanishing field). An immedi- 
ate consequence of (fTO)) is that the two potentials 
take their absolute minimum at the same field value 
P~wp,o = Vo- If follows that the barrier heights 

«(0)-!)(( ( 9o)=i'wp(0)-wwp(^o) (11) 

are identical. Further similarities and differences 
between (fTJ) - © have been studied in [l3l.[Ti|. 



III. NUMERICAL METHODS 

In this section, we summarise the numerical 
techniques used to solve (H} - © and similar to 
high accuracy at a fixed point, and to deduce 
the leading and subleading scaling exponents. 
These include tools for local polynomial expansions 
with various degrees of sophistication (a) - (c), 
tools for initial- value problems (d), and tools for 
solving (two-point) boundary- value problems (e). 
Analytical methods for solving (p} can be found in 
[l6l [TtJ , and are not further elaborated here. 



A. Local behaviour 

For small fields, we use polynomial expansions 
of the effective potential uM . Their convergence 
has been addressed in 

[ISIS HI- We implement 
the expansion in different ways (a) - (c). At a 
fixed point, the potential is ^-symmetric under 
if — ► — (p. Therefore, we expand the effective 
potential in p rather than in (p. 

(a) We Taylor-expand the potential u(p) in field 
monomials about vanishing field ~ p n up to the 
maximum order m in the truncation, 



m ^ 

A ^ Til 



(12) 



Inserting (|12[) into (p} and solving for dtu' = 
leads to coupled ordinary differential equations 
d t X n — Pn({K}) for the couplings Xi, which are 
solved numerically. We note that the m differential 
equations dtX n depend, in general, on m + 2 
couplings up to order X m +2- Therefore, we have to 
provide boundary conditions for the couplings A m +i 
and A m +2, as their values are undetermined by the 
truncated flow. For (fP2)) . we employ A m +i = 
as a boundary condition (in an expansion about 
vanishing field, the dependence on A m +2 drops out), 
which shows good convergence with increasing m; 
see Tab. [1] for numerical results for the couplings at 
the fixed point of (fT|). 

(b) The convergence is further enhanced by ex- 
panding u(p) about non- vanishing field 16, 19]. In 
this case, we write 



u(p) 



m 1 

J2- l x n ( P -x 1 y 



(13) 



where the expansion point p — X\ is defined implic- 
itly by u'(Xi) = c with c a free parameter. From 
(fTJ) we conclude that c can take values between — 1 
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Figure 1: Fixed point potential u(p) from ([l}, ([2]) of 
the Ising universality class for all fields p £ [—00,00]. 
The absolute minimum is at p = 1.814 898 403 687 • • • , 
where the potential is normalised to u(po) = 0. The axes 
are rescaled asi^ 2+\x\ wrt h x ~ p/Po, an d u — > • 
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Figure 2: First derivative of the fixed point potential 
u'(p) from JTJ), 1(2} in the Ising universality class. Nu- 
merically, it reads u' (0) = -0.186 064 249 470 ■■ • at van- 
ishing argument, u' is a monotonously increasing func- 
tion of the field. Same rescaling as in Fig. [T] (see text). 



and 00. This is confirmed by the explicit result, see 
Fig. [21 A natural choice is c = 0, in which case Ai 
denotes the unique potential minimum, u'(Ai) = 0. 
The numerical convergence of the expansion is best 
for small c. In a given truncation to order to, the 
boundary condition X m +i = = X m +2 is used. 
The expansion displays very good convergence 
properties with increasing to [la ]: see Tab. [TJ for 
numerical results for the couplings at the fixed 
point of H|). 

(c) The simple boundary conditions X m +i = 
or A m+ i = = A m+2 used in (a) or (b), respec- 
tively, fail to provide convergent solutions for the 
Wilson-Polchinski flow |5]). The reason behind 
this is that the Wilson-Polchinski flow at small 
fields is more strongly sensitive to the couplings 
at large fields [13j. Identifying a stable solution is 
then more demanding and boils down to providing 
better-adapted boundary conditions for those 
higher order couplings, which are undetermined 
in a given truncation. Equations which fix the 
higher order couplings can be derived from the 
asymptotic behaviour of the scaling solution. An 
algebraic procedure which determines the highest 
order couplings from the asymptotic behaviour of 
@, i.e. from u" — ► for p — ► 00 [and similarly for 
©] has been given in [21(. We have adopted this 
procedure for © (for technical details, see [22|). As 



a result, we find that the polynomial expansion (|12p 
for the Wilson-Polchinski flow ((6]) with appropriate 
boundary condition converges nicely; see Tab. [2] 
for numerical results for the couplings at the fixed 
point of (|5|). Although it is more demanding to 
stabilise the small field expansion based on ([5|) as 
opposed to we note that the Wilson-Polchinski 
potential for small fields differs only mildy from the 
optimised fixed point potential (see Fig. [5] below). 

Polynomial expansions like (fH?)) and (|13p have, 
in general, a finite radius of convergence in field 
space, restricted to field values p m in < p < Pmax- 
Here, we find that p min < and p max > pa- 
in, either of these cases (a) - (c), the flows ([T]) 
and ([5|) transform into a set of m coupled ordinary 
differential equations dtX n = /3«({Aj}) for the co- 
efficient functions A„. The fixed point equations 
(3n({K}) — can be solved to very high accu- 
racy with standard methods. The scaling expo- 
nent v and subleading corrections-to-scaling expo- 
nents are deduced from the eigenvalues of the sta- 
bility matrix M at criticality, My = d^P^g^i—Q. 
Asymmetric corrections-to-scaling, i.e. eigenpertur- 
bations not symmetric under ip — > —ip, can be ob- 
tained as well [23| . 
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coupling 


P = 


u'(p = Xi) =0 


Ai 


-0.186 064 249 470 


1.814 898 403 687 


A 2 


0.082 177 310 824 


0.126 164 421218 


A 3 


0.018 980 341099 


0.029 814 964 767 


A 4 


0.005 252 082 509 


0.006 262 816 384 


As 


0.001 103 954106 


-0.000 275 905 516 



Table 1: The first five couplings at the fixed point of the 
optimised flow {TJ at vanishing field (first column), and 
at the potential minimum (second column). 



coupling 


P = 


w'(p = Ai)=0 


a wp 


-0.228 598 202 437 


1.814 898 403 687 


\WP 


0.187 236 893 730 


0.086 535 420 434 


\WP 


-0.105 930 606 484 


-0.028 253169 622 


Af p 


0.101611201027 


0.015 928 269 983 


\WP 


-0.135 786 295 049 


-0.012 666 298 430 



Table 2: The first five couplings at the fixed point of the 
Wilson-Polchinski flow ([5]) at vanishing field (first col- 
umn), and at the potential minimum (second column). 
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Figure 3: Mass function u + 2p u" of the fixed point po- 
tential (0 in three dimensions, Ising universality class. 
Same rescaling as in Fig. [T] (see text). For p = |y> 2 > 0, 
the mass function is equivalent to v" (tp) from @. 



B. Global behaviour 

Next we discuss two methods (d) and (e) which 
integrate @ or © directly, without using polyno- 
mial approximations. 

(d) Initial value problem. The differential equa- 
tions @ or © are studied as initial value prob- 
lems with boundary conditions given at some start- 
ing point p\. The necessary boundary conditions 
can be obtained within the polynomial expansion 
(a) and (b) in their domain of validity, see Tab. Q] 
for the couplings at vanishing field and at the po- 
tential minimum. The boundary conditions are 

u'(0) = A ls u"(Q) = A 2 =-|A 1 (l + A 1 ) 2 
u'(0) = Af p , u"(0) = Xf p = -| A fP(i + AfP) 

at pi = for the optimised RG flow [l6| and the 
Wilson-Polchinski flow, respectively. Identifying 
the well-defined fixed point solution u'{p) which 
extends over all fields p requires a high degree 
of fine-tuning in the boundary condition Ai [20(. 
Integrating ^ towards larger fields, starting at 
some non- vanishing field |p| ~ 0(1) with boundary 
conditions from (a) or (b) is numerically more 
stable. Following this strategy, we have computed 
in Figs. [TJ [5] and [3] the fixed point potential u, its 
first derivative u' , and the field-dependent mass 
term u' + 2pit", respectively, for all p £ [—00,00]. 



We note that u displays only one global minimum. 
Both u' and u' + 2pu" are monotonously increasing 
functions of the field variable. The numerical 
solution fits the expected analytical behaviour for 
asymptotically large fields. 

(e) Finally, we discuss solvers for two-point 
boundary value problems. The shooting method 
[23 | requires boundary conditions at two points pi 
and P2 in field space. Initial conditions at p\ are var- 
ied until the boundary condition at pi is matched. 
The procedure is iterated until the desired accuracy 
in the solution is achieved 24]. If one is not con- 
strained by the Zi symmetry tp — ► —p and because 
of the potentially singular behaviour of the right- 
hand sides of ([5]) and (O for p — > 0, it is prefer- 
able to implement (U|) and © . Then the differential 
equations at the origin tp = are better under con- 
trol and one may shoot from the origin tp± = (with 
the required symmetry conditions imposed) to large 
fields tp2 ~ Abound j where the asymptotic large-field 
behaviour is imposed. Shooting in the reverse di- 
rection may provide a better accuracy. The large- 
field behaviour of , © and (0J , ([5]) h as p reviously 
been determined in the literature, e.g. [131 ]. The RG 
eigenvalues are deduced from dJ) or §5$ in the vicin- 
ity of the fixed point solution. The implementation 
of (O, ([3]) and (0, © in terms of tp allows a direct 
computation of asymmetric correction-to-scaling ex- 
ponents. 



■5 



iV a , 



rn 



m 



Figure 4: Accuracy of the fixed point, as defined in (I14|l . 
and its dependence on the order of the expansion to; 
from fTJ using method (b) with c = 0. 



Figure 5: Accuracy of the scaling exponent v, as de- 
fined by (115[) . and its dependence on the order of the 
expansion m; from {1} using method (b) with c = 0. 



IV. ERROR CONTROL 

Within the numerical approaches (a) - (e), there 
are several sources for numerical errors, the control 
of which is discussed in this section. 

(i) Within the polynomial expansions (a) - (c), 
the fixed point is determined by seeking simultane- 
ous zeros of all /?- functions. Solutions for d t X n = 
are found to very high precision. To ensure that 
the fixed point is numerically a good approximation 
to the full fixed point solution, we have computed 
\dtu'(p)\ for all p within the domain of validity of the 
polynomial approximation. This serves as a mea- 
sure for the quality of a polynomial expansion to 
order m. We define the accuracy A^cc of the fixed 
point solution to order m as 
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-JV a 



max 

pG[0,p ma 



\dtu'{p)\ 



(14) 



Here, [0, p ma x] denotes a compact neighbourhood of 
the expansion point Ai, which needs not to coincide 
with the potential minimum po . We take p max > Po 
from u'(p max ) + 2/9 max u"(/9 max ) = 1. This notion 
of accuracy is a good measure for how well a local 
Taylor expansion to order m approximates the full 
fixed point potential on [0,p max ]. Using (b) with 
c = 0, we find that (fl"4")) achieves its extremum 
typically at p = 0. In Fig. QJ we display N acc as a 
function of m. Full red (open blue) dots indicate 
that m&xdtu' is positive (negative). The slope is 
approximately 1/2, indicating that an increase in 
the truncation by Am ps 2 increases the accuracy 
in the fixed point solution roughly by a decimal 



place. This serves as an indicator for the reliability 
in the scaling exponents [see (ii)]. 

(ii) Within the polynomial expansions (a) - (c), 
we study the numerical convergence of both the 
fixed point values and the scaling exponents with 
increasing degree of truncation m (l6j . In analogy 
to IT4")) , we define the number of significant figures 
Ax in a fixed point coupling or a critical index X m 
at order m in the expansion as 
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-N x _ 



1 - 



(15) 



where X denotes the full (asymptotic) result. For 
the leading scaling exponent X = v, we display 
N v as a function of m in Fig. O based on the 
expansion (fl~3|) with c = 0. An open blue (full red) 
dot indicates that v(m) is smaller (larger) than 
the asymptotic value. The expansion converges 

roughly in the pattern + H . Note that the 

slope, in comparison with Fig. 21 slightly decreases 
towards larger m. This part of the analysis is 
conveniently performed with Mathematica. The 
accuracy of the matrix inversion (leading to the 
scaling exponents) has been checked independently 
with standard routines from Matlab. 

(iii) Within the polynomial expansions (a) and 
(b), we study the numerical stability of the result 
by varying the expansion point. The radius of 
convergence of polynomial expansions depends on 
the latter. This check serves as an indicator for a 
possible break-down of the expansion at the highest 
orders. We have confirmed that only the rate of 
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Figure 6: Comparison of the fixed point potentials 
v{tp) and their first derivatives v'(<p) from the optimised 
RG (full lines) and the Wilson-Polchinski RG (dashed 
lines). Both potentials are normalised to v(ipo) = 0, 
where they take their absolute minimum. In either case, 
ipa = y/2p^ = 1.905 202 563 344 • • ■ . 



Figure 7: Comparison of the fixed point potentials 
v(tp) from the optimised RG (full line) and the Wilson- 
Polchinski RG (dashed line) for all fields tp £ [0, oo]. As 
in Fig. [6] the potentials are normalised to v(ipo) = 0. 
For display purposes, the axes are rescaled as x — > 2 +\x\ 
with x — tp/tpo, and v — > r^—. 



convergence depends on the expansion point, but 
not the asymptotic result. 

(iv) The numerical approach (d) is checked in 
several ways. First, starting at intermediate fields, 
the numerical accuracy in the integration can be 
made large. Second, the domain of validity for (d) 
and (a) - (c) overlap. This allows for a quantitative 
cross-check. Third, for large fields \p\ — * oo, the 
results from (d) match the expected asymptotic 
behaviour. Fourth, we have cross-checked the result 
with a two-point boundary value routine. 

(v) The shooting method (e) is controlled in sev- 
eral ways. First, the numerical precision for the so- 
lution on [0, pbound] is only limited by the machine 
precision. We use standard routines under Fortran 
f 77 with double precision. The boundary condition 
at large fields has not to be known to very high ac- 
curacy to achieve a reliable result. Stability in the 
result is confirmed by varying /Obound as well as the 
boundary condition at pbound- These procedures are 
applied for both ([2]) and ([6]). The shooting method is 
checked independently using a different solver under 
Maple, based on a different boundary value inte- 
gration algorithm. 



V. RESULTS 

Our results for the scaling exponents from ([I]) 
and ([7]) are given in Tab. [3] and |U All digits are 
significant, except for possible rounding effects. 
Earlier findings with a lower level of accuracy 
[IS IM H are full y consistent with Tab. O Specif- 
ically, in Tab. [3l the exponents from the optimised 
RG flows have been obtained using (a), (b), and 
(e), whereas those for the Wilson-Polchinski flow 
follow from (c) and (e). The method (b) has 
been used up to the order to = 32, leading to 
roughly 11 significant digits in the result for v, and 
roughly one significant digit less with increasing 
order for the subleading exponents, see [IS HH- 
Method (e) has been pushed to an accuracy of 14 
digits for all scaling exponents. The asymmetric 
corrections-to-scaling exponents in Tab. 3] have 
been computed using (e), again to an accuracy of 
14 digits. Again, all earlier results to lower order in 
the accuracy (23|, are consistent with Tab. [U 

In Figs. Q] - 02 we have plotted our results for 
the fixed point potential u(p) from ([2]) and its 
derivatives for all p G [—00,00]. It is noteworthy 
that the solution extends to all negative p, which 
in terms of the field ip corresponds to the analytical 
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exponent 


optimised RG 


Wilson-Polchinski RG 


V 


0.649 561 773 880 65 


0.649 561 773 880 65 


w 


0.655 745 939 193 3 


0.655 745 939193 35 


UJ2 


3.180 006 512 059 2 


3.180 006 512 059 2 


U)3 


5.912 230 612 747 7 


5.912 230 612 747 7 


U>4 


8.796 092 825 414 


8.796 092 825 414 


^5 


11.798 087 658 337 


11.798 087 658 336 9 


LOd 


14.896 053175 688 


14.896 053175 688 



Table 3: The leading scaling exponent v and the first six 
subleading eigenvalues for the Ising universality class in 
three dimensions. 



continuation to purely imaginary field values. For 
p — > — oo, the solution approaches the convexity 
bound where v! + 2pu" — ► — 1. In the flow equation 
this limit is potentially singular as it corre- 
sponds to a pole. However, the explicit solution 
for d t u' — shows that the pole is supressed, 
leading to a well-behaved solution even in this 
limit. In the physical domain where p > 0, the 
fixed point solution stays far away from the pole, 
since min p >o(u' + 2pu") = u'(0) > — 1. General 
(non-fixed point) solutions to |T]) approach the 
pole only in a phase with spontaneous symmetry 
breaking where V" fc "(</>) is negative. Then the pole 
ensures convexity of the physical potential Vk=o(<f>) 
in the infrared limit where V£' =0 ((j>) > [271 ] - 

In Figs. [5] and [7] we compare the fixed point 
potentials v(ip) from the optimised RG and the 
Wilson-Polchinski RG for small and large fields. 
We have checked that u(p) in Fig. [1] matches 
v(<p) in Fig. as it should. Both potentials have 
their absolute minimum at the same field value 
P™ P = Po = 1.814 898 403 687; see Tab. ED and H 
It is also confirmed that the potentials display the 
same barrier height v(0) — v(tpo) as expected from 
(fTTj) . For large fields, the potentials scale differently 
with the fields fTl. 



We have checked numerically that the fixed point 
solutions to J2J) and ([6]) are related by a Legendre 
transform, see ©, ([1(3) . Consequently, at vanishing 
field, the couplings from Tab. Q] are related to the 
Wilson-Polchinski couplings in Tab. [2] by 



Ai = 



WP 



A I 



WP 



1 _ a wp - 

Ai 
1 + Ai ' 



A, 



\ W P 

2 — 



(1 - AfP) 4 
A 2 



(1 + XxY 



(16) 
(17) 



Similar expressions are found for the higher order 
couplings. The numerical results in Tab. [1] and [2] 
confirm these relations within the present accuracy. 



exponent 


optimised RG 


Wilson-Polchinski 


Q 


1.886 703 838 0914 


1.886 703 838 0914 


UJ2 


4.524 390 733 670 8 


4.524 390 733 670 8 


UJ 3 


7.337 650 643 354 


7.337650 643 354 4 


U>4, 


10.283 900 724 026 


10.283 900 724 026 



Table 4: The first four asymmetric correction-to-scaling 
exponents for the Ising universality class in three dimen- 
sions. 



Our value \f p = -0.228 598 202 43702 for the 
dimensionless mass term squared at vanishing field 
deviates at the order 10~ 5 from the corresponding 
value -0.228 601 293 102 given in [H]. 

In summary, we have confirmed the equivalence 
of an optimised RG and the Wilson-Polchinski RG 
to the order 10 -14 in the universal indices, and in 
the Legendre-transformed fixed point solutions. 



VI. COMPARISON 

Next we compare our results with those based 
on Dyson's hierarchical model 0. To leading order 
in the derivative expansion, it has been proven by 
Felder [ll[ that the hierarchical RG is equivalent 
to the Wilson-Polchinski RG, and hence to the 
optimised RG. Then scaling exponents should come 
out identical. With high-accuracy numerical data 
at hand, we can test this assertion quantitatively. 

In hierarchical models, the block-spin trans- 
formations are characterised by a decimation 
parameter I > 1, where I — > 1 refers to Felder 's 
limit of continuous transformations. Numerical 
studies with £ = 2 1 / 3 « 1.26 have been performed 
in [H |U EH H HI El, El- In Tab. E we 
compare our results with previous computations 
of v from the Wilson-Polchinski RG [H HH, the 
optimised RG 113. Rol . and the hierarchical RG 
pj, ES EH ElUElS El- The most advanced 
hierarchical RG (HRG) and functional RG (FRG) 
studies agree amongst themselves at least to the 
order 10 -6 and 10 -12 , respectively. We emphasize 
that the functional RG results clearly deviate 
from the hierarchical RG results, although the 
difference (^frg - ^hrgV^frg ~ -1-3 X 10~ 5 is 
quantitatively small. 

In Tab. [S] we compare the best results from 
the functional RG (i.e. both the Wilson- 
Polchinski RG and the optimised RG) with 
those from the hierarchical RG for various other 
indices, some of which have been obtained 
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Wilson — Polchinski optimised RG hierarchical RG 


exponent 


functional RG 


hierarchical RG 


0.649 c 0. 64956 s 0.6496" 


V 


0.649 561773 880 65 


0.649 570 e ' Li 


0.6496 d 0.649562' 1 0.64957 ! ' 


LU 


0.655 745 939 193 3 


0.655 736 i 


0.649561773881 0.649561773881 0.649570 e,/ ' i 


7 


1.299 123 547 7613 


1.299 140 730 159 e 




A 


0.425 947 495 477 4 


0.425 946 859 881 e 


Table 5: Comparison of the scaling exponent v for the 


a 


0.051314 678 358 05 


0.051 289 i 



Ising universality class in three dimensions (see text) 
Data from this work, and from o) [301, 6) [3jJ , c) [28|] 



d) [28 



/) 33], g) 12 



Table 6: Comparison of universal indices v, ui, 7 = 2 v, 
A = to v and a = 2 — ?>v (for 77 = 0) from the func- 
tional RG (this work) and the hierarchical RG (same 
referencing as in Tab. [5}. 



with up to 12 significant figures [32[ (see also 
[H|). The variations (X FRG - ^hrgV^frg 
in the exponents X = v, oj,7, A, a read 
(-1.3,1.5,-1.3,0.15,5.0) x 10~ 5 , respectively. 
Hence, we confirm that the findings disagree, 
beginning at the order 10~ 5 . Given Felder's proof 
of equivalence for the limit t — > 1 [ll[ , we conclude 
that the hierarchical model diplays a dependence 
on the decimation parameter £, with a tiny slope of 
the order of 10 -6 for variations in the exponents. 



VII. DISCUSSION AND CONCLUSIONS 

Equivalences between non-linear differential or 
difference equations, and more generally functional 
relationships between different implementations of 
the renormalisation group, often allow for a deeper 
understanding of both the underlying physics and 
the adequacy and efficiency of the methods at 
hand. It has been proven previously that three dif- 
ferent implementations of Wilson's renormalisation 
group - the Wilson-Polchinski RG, an optimised 
version of the effective average action RG, and 
Dyson's hierarchical RG in the limit of continuous 
block-spin transformations - are equivalent in the 
local potential approximation. This implies that 
the corresponding non-linear RG equations carry 
identical universal content, e.g. identical scaling 
exponents. 

Quantitatively, these equivalences have been 
confirmed up to the order 10~ 4 in the literature, 
but a discrepancy at the order 10~ 5 has recently 
been pointed out. This last observation relies on 
studies within an optimised RG and Dyson's hier- 
archical RG with decimation parameter I — 2 1 / 3 , 
where results with a sufficiently high accuracy 
are available, while previous results from the 
Wilson-Polchinski RG were in agreement with 
either of them. Here, we have closed this gap in the 



literature by computing scaling exponents for the 
Ising universality class from the Wilson-Polchinski 
RG and the optimised RG with an unprecedented 
accuracy of the order 10~ 14 . We confirm their 
equivalence, first conjectured in fl2 ], based on the 
leading, subleading and asymmetric correction-to- 
scaling exponents. Equally important, our central 
numerical results are obtained in several ways, 
and furthermore independently from both the 
Wilson-Polchinski and the optimised RG flow. We 
conclude that their equivalence is rock solid for all 
technical purposes, and in full agreement with the 
explicit proof given by Morris [14| . 

In contrast to this, we now have clear indications 
for a non-equivalence of our results with those from 
Dyson's hierarchical model for discrete block-spin 
transformations. In the 3d Ising universality 
class, all scaling exponents from the functional 
RG, i.e. from both the Wilson-Polchinski and the 
optimised RG, differ systematically at the order 
10~ 5 from high accuracy studies based on Dyson's 
hierarchical RG. Given the high degree of accuracy 
in the critical indices from the present and earlier 
studies, and the fact that the scaling exponents 
in all three approaches have been obtained from 
several independent numerical implementations 
and collaborations, it is unlikely that the dis- 
crepancy originates from numerical insufficiencies. 
Furthermore, these differences are absent in the 
limit of continuous block-spin transformations 
[ill ]. We conclude that the tacit assumption of 
^-independence in t he scaling exponents of the 
hierarchical model 



15] 



cannot be maintained. 
Rather, the hierarchical model carries an inherent 
^-dependence, similar to a dependence previously 
observed in [3(| (see also [53] )• In functional RG 
approaches, scheme-dependences related to the 
underlying Wilsonian momentum cutoff can arise 
in truncations [2j. Here, powerful control and opti- 
misation mechanisms are available which decrease 
truncational artefacts and increase convergence 
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towards the physical theory 0, H, [38j]. It will be 
interesting to see if the ^-dependence of hierarchical 
models or its extensions can be understood along 
similar lines. 

For our numerical work, we have developed meth- 
ods to solve non-linear eigenvalue problems with 
high precision. The combined use of local polyno- 
mial expansions techniques with global integration 
methods allows for an efficient determination 
of the full fixed point solutions with reasonable 
numerical efforts. We have also put some emphasis 



on a reliable error control, in view of the high 
accuracy aimed for in the universal eigenvalues. 
We expect that this combination of techniques will 
be equally useful for other non-linear problems in 
mathematical physics, including e.g. the derivative 
expansion to higher order. 
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